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Topology, fermionic zero modes and flavor singlet correlators in finite 
temperature QCD. 
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We extend our earlier results concerning the breaking of the £/a(1) symmetry in finite temperature QCD. In 
particular, we use topological argments to investigate the chiral limit. 



Among the non-perturbative properties of 
quarks, those associated with the topology of the 
gauge fields occupy a special place. This fol- 
lows from the Atiyah-Singer index theorem which 
guarantees the existence of fermionic zero modes 
on configurations with non-trivial topology. As 
a result, non-perturbative issues related to the 
axial anomaly can generally be studied with en- 
hanced control. In the context of finite temper- 
ature QCD, one such issue relates to the realiza- 
tion of the anomalous U(l) axial symmetry. If 
this symmetry were to be effectively restored at 
the critical temperature, the increased number of 
soft degrees of freedom would modify the nature 
of the phase transition |lj]. Needless to say, only 
the lattice can really hope to tackle this question, 
at this time. Unfortunately, the lattice formu- 
lation also introduces artifacts and there is for 
example no exact index theorem in this context. 
Instead, what is observed with the actions cur- 
rently used is the "zero mode shift" phenomenon 
|2| (i.e. modes with small but non-zero eigen- 
value and which are not exactly chiral) . What we 
would like to show here is that even in this case, 
enough is recognizable of the continuum behavior 
that meaningful conclusions can be drawn. We 
find that the Ua(1) symmetry is not restored at 
the phase transition but only at somewhat higher 
temperatures. At the same time, it is clear that 
detailed quantitative questions could only be an- 
swered with the use of an improved action. 
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Some of the results presented here are an ex- 
tension of preliminary data introduced last year 
||. Many of the trends are confirmed. In addi- 
tion, we discuss the issue of reaching the chiral 
limit. We have measured the screening correla- 
tors in the scalar and pseudoscalar flavor singlet 
and flavor triplet sectors. From there, one can 
investigate the symmetry restorations by look- 
ing for degeneracies, either in the fitted screening 
masses or in the susceptibilities (integrated corre- 
lators). These computations were carried out on 
a set of configurations from the HTMCGC collab- 
oration These were generated on a lattice of 
size 16 3 x 8 with 2 flavors of dynamical staggered 
fermions with mass ma — 0.00625. In addition 
to those, we have generated a new set of config- 
urations at [3 — 5.4875 which is just above the 
crossover temperature at this value of the quark 
mass. Our results for the screening masses are 
presented in figure 1. The new measurements 
at j3 — 5.4875 confirm the trend observed ear- 
lier, namely that the cr becomes light close to the 
transition while the 8 does not. It is worth men- 
tioning here that some of the fits presented in 
fig.l may have large systematic errors (only the 
statistical errors are included in the figure) . This 
is due in part to the small extent of the lattice 
and, in the case of the cr, to the additional diffi- 
culties associated with the measurements of dis- 
connected quark loop correlators. Nevertheless, 
it seems to us that the trends observed in fig. 1 
are rather clear. In addition, they are consistent 
with our measurements of the susceptibilities (as 
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Figure 1. Screening masses in two flavor QCD 



well as those by other groups j5j), in that the peak 
of the scalar susceptibility is seen in the discon- 
nected and not the connected contribution (for 
summary plots of earlier results see ref ||). If 
this situtation persists in the chiral limit, it will 
imply a non-restoration of the E/a(1) symmetry 
in the high temperature phase of QCD. Namely, 
the 7r and a become light close to the transition, 
but not the S. 

In the case of the susceptibilities, extrapola- 
tions to the chiral limit at fixed (3 have been at- 
tempted by two groups but the results appear 
ambiguous For example, it is seen that lin- 

ear fits generally have a better \ 2 than fits to a 
quadratic curve which would be expected on the- 
oretical grounds. Here we would like to argue 
that these difficulties are due to lattice artifacts 
and will propose a method which allows us to ex- 
tract information about the chiral limit without 
carrying out an explicit extrapolation. The idea 
is to identify the reason for the non-vanishing of 
the disconnected susceptibility in the chiral limit. 
This is done by making use of the spectral decom- 
position of the quark propagator and the resulting 
formulas: 
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where til (tir) are respectively the number of 
left (right) zero- modes of the Dirac operator. In 
terms of these the scalar and pseudoscalar dis- 
connected susceptibilities are written x dls = [< 
Q 2 > -(< Q >) 2 }/V and X f s =< Ql > /V. 
Now it is clear that only the contribution from 
fermionic zero-modes will survive at m = and 
that the chiral limit will be smooth. Remember 
that for Nf = 2 the fermionic determinant picks 
up a factor of m 2 for each zero mode. On config- 
urations with topological charge one, the factors 
of m will cancel from measurements of xfe and 
leave a well defined non-zero result. This is what 
we expect in the continuum, on the lattice how- 
ever things will be more complicated. First, there 
are no exact zero modes ( except on a subspace 
of measure zero ). Then: 



\Latt 



2m 



Q 



Latt 



— Ea>0 X 2 +m 2 
J2 2m<-0 A |r 5 |^A> 



;A>0 ' 



(4) 
(5) 



S(x,y) 



E 



A i\-\-m 



(i) 



where we have used the symmetries of the stag- 
gered action and T§ is the (four-link) lattice 75 
operator. It is then clear that on the lattice, dis- 
connected susceptibilities will vanish in the chiral 
limit and be proportional to m 2 for small m. We 
believe that it is this lattice artifact which makes 
a chiral extrapolation from the lattice data ex- 
tremely difficult UU . 

We therefore propose a different method for 
studying the m — > limit. The idea is to iden- 
tify "continuum behavior" in the lattice data (and 
then use what wc know about the continuum the- 
ory). Together with the Atiyah-Singer index the- 
orem, (3) implies mQ^ = Qt op - We would like to 
see the same connection with topology on the lat- 
tice. The formula (3) comes about because in the 
continuum < ^a^sI^a > is either ±1 if A = or 
otherwise. The lattice expression (5) will closely 
match the continuum if we see on the lattice a 
clear correlation between small eigenvalues and 
large r\ =< "0a |Ts |-0a > (i-e. if we can identify 
chiral modes). Figure 2 is a plot of r\ versus A 
for the lowest 8 positive eigenvalues on each con- 
figuration of our sample at (3 — 5.55. Many of 



3 



Np-2 6-5.55 m n a-0. 00625 



0.150 



0.125 



0.100 - ° 



Crl 0.075 - n m □ ° 



050 - U ° *nS° D i 
0.025 h fn a n n > 



0.000 




■1 0.02 0.03 

X 

Figure 2. (r, A) diagram at (5 = 5.55 
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Figure 3. disconnected 75 susceptibility 



the small eigenvalues are indeed associated with 
a larger value of r\. Of course, the quality of 
the correlation depends on how strong the cou- 
pling is ( It is for example not as clear here as 
it was in quenched QCD at [3 = 6.2 @ ). The 
second step in the argument consists in showing 
that these low modes (the would be chiral modes 
of the continuum theory) indeed saturate the lat- 
tice results. In figure 3 we compare the full result 
for xf ls with the one obtained by truncating (5) 
to the lowest 8 modes. The agreement is per- 
fect at (3 = 5.55 and deteriorates somewhat in 
the rest of the symmetric phase (in part because 
of lattice artifacts). Putting everything together, 
we find that xt ls ^ s non-zero in the symmetric 
phase and that the topological origin of this re- 
sult guarantees that it will remain non-zero in 
the continuum chiral limit. We therefore con- 
clude that the Ua(1) symmetry is not restored 
at the phase transition ^. Our results also in- 
dicate the importance of investigating improved 
lattice actions which better satisfy the index the- 
orem. These would allow the continuum behav- 
ior to appear even more clearly than here. Also, 



3 Of course, this conclusion only concerns the volume stud- 
ied so far, namely V = (2/T) 3 . The actual volume depen- 
dance remains to be investigated. 



lower values of the quark mass would be allowed 
before the zero-mode shift phenomenon becomes 
a problem. Certainly, such actions will be nec- 
essary to obtain quantitative results. In fact, 
some of the formulations that one would like to 
consider have been mentioned at this conference 
(domain wall fermions, perfect action, ...). Be- 
yond that, there are many other questions which 
can be investigated with the knowledge of the low 
lying eigenvalues and eigenvectors. Correlations 
between the eigenvalues is one of them. The re- 
lation with instantons and the possible existence 
of instanton- anti-instanton molecules is another. 
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